Stable Fermion Bag Solitons in the Massive Gross-Neveu Model: Inverse Scattering 

Analysis 
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Formation of fermion bag solitons is an important paradigm in the theory of hadron structure. We 
study this phenomenon non-perturbatively in the 1 + 1 dimensional Massive Gross-Neveu model, 
in the large N limit. We find, applying inverse scattering techniques, that the extremal static 
bag configurations are reflectionless, as in the massless Gross-Neveu model. This adds to existing 
results of variational calculations, which used reflectionless bag profiles as trial configurations. Only 
reflectionless trial configurations which support a single pair of charge-conjugate bound states of 
the associated Dirac equation were used in those calculations, whereas the results in the present 
paper hold for bag configurations which support an arbitrary number of such pairs. We compute 
the masses of these multi-bound state solitons, and prove that only bag configurations which bear 
a single pair of bound states are stable. Each one of these configurations gives rise to an 0(2N) 
antisymmetric tensor multiplet of soliton states, as in the massless Gross-Neveu model. 

PACS numbers: ll.10.Lm, 11.15.Pg, ll.10.Kk, 71.27.+a 



INTRODUCTION 

An important dynamical mechanism, by which funda- 
mental particles acquire masses, is through interactions 
with vacuum condensates. Thus, a massive particle may 
carve out around itself a spherical region .1] or a shell 
in which the condensate is suppressed, thus reduc- 
ing the effective mass of the particle at the expense of 
volume and gradient energy associated with the conden- 
sate. This picture has interesting phenomenological con- 
sequences pit!- 

This dynamical distortion of the homogeneous vacuum 
condensate configuration, namely, formation of fermion 
bag solitons, was demonstrated explicitly by Dashen, 
Hasslacher and Neveu (DHN) many years ago, in their 
study of semiclassical bound states in the 1 + 1 dimen- 
sional Gross-Neveu (GN) model [5j- 

Fermion bags in the GN model were discussed in the 
literature several other times since the work of DHN, us- 
ing alternative methods 0, 0, @ . For a review on these 
and related matters (with an emphasis on the relativis- 
tic Hartree-Fock approximation) see Q- For a more re- 
cent review of static fermion bags in the GN model (with 
an emphasis on reflectionless backgrounds and supersym- 
metric quantum mechanics) see [lfj. The large-iV semi- 
classical DHN spectrum of these fermion bags turns out 
to be essentially correct also for finite N, as analysis of 
the exact factorizable S-matrix of the GN model reveals 

mi. 

A variational calculation of these effects in the 1 + 1 
dimensional massive generalization of the Gross-Neveu 
model, which we will refer to as MGN, was carried in 
[l2| a few years ago, and more recently in [l^. In this 
paper we study static fermion ba gs i n the MGN model 
using inverse-scattering formalism |14j , thus avoiding the 



need to choose a trial variational field configuration. The 
present work is thus a natural extension of the inverse- 
scattering analysis carried out by DHN to the massive 
case. 

The MGN model is defined by the action 
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describing N self interacting massive Dirac fermions i\) a 
carrying a flavor index a = 1, . . . , JV, which we promptly 
suppress. 

An obvious symmetry of with its TV Dirac spinors is 
U(N). Actually, 0) is symmetric under the larger group 
0(2N) (see also Section 1 of lOj). The fact that the 
symmetry group of is 0(2N) rather than U(N), in- 
dicates that it is invariant against charge-conjugation. It 
is easy to see this in a concrete representation for 7 ma- 
trices, which we choose as the Majorana representation 

00 
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(Henceforth, in this paper we will use this representation 
for 7 matrices in all explicit calculations.) 

In the representation (J2J), charge-conjugation is real- 
ized simply by complex conjugation of the spinor 



^(x) = ij)*{x). 



(3) 



Thus, if ip 
tion 



t u(x) is an eigenstate of the Dirac equa- 



[i$ - cr(x)] V = , 



(4) 
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with energy ui (assuming time independent cr(x)), then 



"(x) is an energy eigenstate of (@J, with 



energy —u. Therefore, the MGN model QJ, like the GN 
model, is invariant against charge conjugation, and en- 
ergy eigenstates of (@J come in ±u> pairs. 

As usual, the theory (Q can be rewritten with the help 
of the scalar flavor singlet auxiliary field <r(x). Also as 
usual, we take the large N limit holding A = Ng 2 fixed. 
Integrating out the fermions, we obtain the bare effective 
action 

S[a] = - ^2 y d2x (°" 2 - 2Mcr ) - iN Tr lo S - °") ■ 

(5) 

Noting that 75 (ip — <r) = — (ip + 17)75, we can rewrite 
the Tr log(i^-tr) as ^Tr log - o){i$ + a)]. In this 

paper we focus on static soliton configurations. If a is 
time independent, the latter expression may be further 
simplified to f J" [Tr log(/i+ - w 2 ) + Tr bg(ft_ - w 2 )] 
where 
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and where T is an infra-red temporal regulator. 

As it turns out, the two Schrodinger operators h± are 
isospectral (see Appendix A of [13 and Section 2 of @) 
and thus we obtain 



S[<t] = J d 2 x (a 2 -2Ma) 

OO 

- iNT [ ^Tr log(/i_ - w 2 ) . (7) 
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In contrast to the standard massless GN model, the 
MGN model studied here is not invariant under the Z2 
symmetry ip — > 75^, <r — > —a, and the physics is cor- 
respondingly quite different. As a result of the Z2 de- 
generacy of its vacuum, the GN model contains a soliton 
(the so called CCGZ kink HHHESEi) in which the 
a field takes on equal and opposite values at x = ±00. 
The stability of this soliton is obviously guaranteed by 
topological considerations. With any non-zero M the 
vacuum value of a is unique and the CCGZ kink be- 
comes infinitely massive and disappears. If any soliton 
exists at all, its stability has to depend on the energetics 
of trapping fermions. 

Let us briefly recall the computation of the unique vac- 
uum of the MGN model. We shall follow [13 • For an 
earlier analysis of the MGN ground state (as well as its 
thremodynamics), see [lij . Setting a to a constant we 
obtain from J7J) the renormalized effective potential (per 
flavor) 
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A( M ) 
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where fj, is a sliding renormalization scale with A(/i) = 
Ng 2 (p) and M(/j) the running couplings. By equating 



the coefficient of a 2 in two versions of V, one defined with 
\x\ and the other with /i2, we find immediately that 



_L l_ = Ilo ^ 

A(/il) A(/i 2 ) TT jl 2 



(9) 



and thus the coupling A is asymptotically free, just as in 
the GN model. Furthermore, by equating the coefficient 
of a in V we see that the ratio is a renormalization 
group invariant. Thus, M and A have the same scale 
dependence. 

Without loss of generality we assume that M(jjl) > 
and thus the absolute minimum of JHJ, namely, the 
vacuum condensate m = (a) , is the unique (and positive) 
solution of the gap equation 



dV 
da 
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Referring to (JTJ, we see that m is the mass of the 
fermion. Using 10, we can re- write the gap equation as 
which shows manifestly that to, an observ- 



M( M ) 



A(m) 

able physical quantity, is a renormalization group invari- 
ant. This equation also implies that M{m) — m, which 
makes sense physically. 

Fermion bags correspond to inhomogeneous solutions 



of the saddle-point equation 
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0. In particular, 



8a(x,t) 

static bags <j(x) are the extremal configurations of the 
energy functional (per flavor) 



£[a(x)] 



SHx)] 
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subjected to the boundary condition that a(x) relaxes to 
its unique vacuum expectation value to at x = ±00. 

The rest of this paper is organized as follows. In the 
next section we express £[cr(x)] in terms of the scattering- 
data. DHN have already done most of the work for us, 
save for one crucial piece, specific to the MGN model: 
an expression for (a(x) — to) dx in terms of the scat- 
tering data. In this paper we derive this relation. The 
details of this derivation are relegated to the Appendix. 
We then complete the task of writing down £[<t(x)] in 
terms of the scattering data. We prove that the static 
extrema of £[c(x}] are reflectionless, as in the mass- 
less GN model |4j, and calculate their masses. Such 
an extremal configuration, considered as a scalar back- 
ground in the Dirac equation, typically supports a num- 
ber K of pairs of charge conjugate bound states at en- 
ergies ±w„ (n — 1,2,..., K) which bind fermions and 
antifcrmions. Each one of these tu n s depends only on the 
total number v n of fermions and antifcrmions it binds, 
thus giving rise to an 0(2N) iank-v n antisymmetric ten- 
sor multiplet of soliton states. As it turns out, the mass of 
such a soliton is a convex function of the v n s. In Section 
3, we invoke the convexity of the soliton's mass formula 
and prove that only solitons which support a single pair 
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of bound states (i.e., K = 1) are stable against decaying 
into lighter solitons. These are precisely the configura- 
tions studied variationally in 0, . 



STATIC CONFIGURATIONS AND 
INVERSE-SCATTERING ANALYSIS 

A fermion bag is essentially a trap for fcrmions. Ev- 
idently, an apriori specification of a static fermion bag 
configuration should indicate the list of bound states 
and enumerate the fermions it traps. To be specific, we 
shall evaluate in this section the energy functional 
of a static configuration o~(x), obeying the appropriate 
boundary conditions at spatial infinity, which supports 
K pairs of bound states of the Dirac equation at energies 
±w„, n = 1,...,K (where, of course, w n < to 2 ). The 
bound states at ±u> n are to be considered together, due 
to the charge conjugation invariance of the GN model. 
Due to Pauli's exclusion principle, we can populate each 
of the bound states ±aj„ with up to N fermions. In such a 
typical multiparticle state, the negative frequency state 
is populated by N — h n fermions and the positive fre- 
quency state contains p n fermions. In the parlance of 
Dirac's hole theory, this is a many fermion state, with 
p n particles and h n holes occupying the pair of charge- 
conjugate bound states at energies ±w„. We shall refer 
to the total number of particles and antiparticles trapped 
in the n-th pair of bound states 



V n = Pn + K 



(12) 



as the valence, or occupation number of that pair. 

From @ and we obtain the bare energy functional 

as 

oc 

£[<j{x)] = J dx [V(x) - 2M(i{x)] 

— QO 

oo 

- / ^Trlog[-d 2 + V(x)-^ 2 ] (13) 



where 



V(x) =a 2 (x) -a'(x). 



(14) 



(Here we used dx o~'(x) = by invoking the bound- 
ary conditions <r(x) to.) The expression (|13f) is di- 
vergent. We regulate it, as usual, by subtracting from it 
the divergent contribution of the vacuum configuration 
a 2 = to 2 and by imposing a UV cutoff Aonw. Thus, the 
regulated bare energy functional associated with a(x) is 



2A 



dx [V(x) - to 2 - 2M(a(x) - m)] 



A 

4 J £ [Trlog^-^-Trlog^-c 2 )] (15) 



where 



h = -d 2 

ll VAC U X 



(16) 



is the hamiltonian corresponding to the vacuum config- 
uration. We are not done yet, as the integrals over to in 
(|15|l still diverge logarithmically with the UV cutoff A. 
However, this logarithmically divergent term is precisely 
the one that should be added to the bare 1/ A in the first 
term in (|15f) in order to obtain the renormalizcd coupling 
appearing in @ HEE3- 

Now that we have written the energy-functional i|15|) of 
a static configuration, or a fermion bag, our next step is 
to identify those fermion bags on which l|15l) is extremal. 

The energy functional (|15f) is, in principle, a compli- 
cated and generally unknown functional of o~(x) and of 
its derivatives. Thus, the extremum condition 



S£[a] _ 



8a{x) 



= 0, 



as a functional equation for o~(x), seems intractable. The 
considerable complexity of the functional equations that 
determine the extremal <j(x) configurations is the source 
of all difficulties that arise in any attempt to solve the 
model under consideration. 

DHN found a way around this difficulty in the case 
of the GN model |^. They have used inverse scatter- 
ing techniques 0] to express the energy functional £ [a] 
H15JI in terms of the so-called "scattering data" associ- 
ated with, e.g., the hamiltonian /i_ in 10 (and thus with 
o~(x)), and then solved the extremum condition on 115|) 
with respect to those data. 

The scattering data associated with are [l| the re- 
flection amplitude r(k) of /i_ at momentum fc, the num- 
ber K of bound states in and their corresponding 
energies < u; n < to 2 ,(n = 1, . . . , K), and also addi- 
tional K parameters {c„}, where c n has to do with the 
normalization of the nth bound state wave function ?/>„ of 
More precisely, the nth bound state wave function, 
with energy w 2 , must decay as ip n (x) ~ const 
as x — > oo, where 



< K r , 



TO" 



(17) 



If we impose that ip n {x) be normalized, this will deter- 
mine the constant coefficient as c„. (With no loss of gen- 
erality, we may take c„ > 0.) Recall that r(— k) = r*(k), 
since the Schrodinger potential V(x) is real. Thus, only 
the values of r(k) for k > enter the scattering data. The 
scattering data are independent variables, which deter- 
mine V(x) uniquely, assuming V(x) belongs to a certain 
class of potentials which fall-off fast enough towards in- 
finity. 

Since the MGN does not bear topological solitons, nei- 
ther h- nor h+ can have a normalizable zero energy 
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eigenstate [l7|. (See also Section A. 1.1 in Appendix A 
of 0.) Thus, all the u n are strictly positive. 

We can apply directly the results of DHN in order to 
write down that part of (|15f) which is common to the 
MGN and GN models, i.e., I|15l) with its term propor- 
tional to M removed, in terms of the scattering data. 
In what follows we briefly summarize their results. (See 
Sections 2, 3, and Appendix B of for details.) 

Using the trace identities of the spectral theory of h~, 
they were able to show that 

oo 

h[r(k),{Kn}] = ~ J dx[V(x)-m 2 } 

— OO 

oo K 

= 2^/ l°g[l-K*0| 2 ] dk+jJ^Kn (18) 

-oo ™ =1 

This takes care of the first term in l|15|) . Onto the spec- 
tral determinants: These encode the contribution of the 
Fermi fields to the energy of the bag relative to the vac- 
uum. To account for it correctly, we should put our field 
theory in a big spatial box of length L so as to make the 
spectrum discrete, which will enable us matching each 
mode of fermion fluctuations around the vacuum, with 
its counterpart, obtained as u{x) is turned on adiabati- 
cally. We shall take the limit L — > oo only in the end. 
Thus, we obtain the Fermi field part of the energy (per 
flavor) as 

A 

£f = % j ^ [Trlog^-^^-Trlog^-c 2 )] 



,VAC\ 



K K 



71=1 



n=l 



(19) 



The first sum runs over all positive-energy scattering 
states, where u> Q [<7(a;)] is the energy of the scattering 
state to which the Fermi mode energy iv^ AC flows to 
as the vacuum configuration a = m is deformed adia- 
batically to <j(x). The second sum in (|19|) accounts for 
the first K scattering states above the threshold u> = to 
which migrate into the gap to become the (positive en- 
ergy) bound states ui n as cr(x) is switched on. In the 
limit L — > oo their energies are indistinguishable from 
to. Note the minus sign in front of these two sums, as 
appropriate for fermion zero-point energy. The last sum 
in (|19JI is evidently the contribution of valence fermions 
and antifermions trapped inside the bag at the bound 
states ±w„. By carefully counting scattering modes in 
the box, DHN arrived at the fairly standard result 



VAC\ 



5{u)duj (20) 



where S(u>) is the scattering phase-shift. Then, chang- 
ing to momentum space (fc 2 + to 2 = ui 2 ), and using a 
dispersion integral representation for 5(k), DHN derived 
that 



h[r(k),{K n }] ee J 5{uj)duj = 

m 

kdk „„ 7 log [1 - \r(q)\ 2 ] 
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kdk 



(fc 2 + to 2 ) 1 / 



K 

arctan— . 



dq 



(21) 



The second integral in 121J1 can be calculated explicitly. 
In the limit A/m >> 1 we obtain 



K r 



I 2 = 2J2 



71=1 



K n log 



/ 2Ae\ 7tto oj n 

h Lu n arctan — 

\ to J 2 K n 



(22) 

DHN's results, Eqs. JIHJ - l|2I}, correspond to all 
terms of ()15(l . save for the term proportional to M, 
J- ( C7 ( x ) — m ) dx. This integral cannot be expressed in 
terms of the scattering data based on the trace identities 
of the Schrodinger operator h— discussed in Appendix B 
of £| . Evidently, new analysis is required to obtain its 
representation in terms of the scattering data. Happily 
enough, we were able to obtain such a representation. In 
the Appendix we prove that 



(<y(x) — to) dx = 
" log[l-|r(g)| 2 ] 
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log 



im — q 



m — Kr, 



TO + K Tl 



dq 



(23) 



Combining (JTSJ- l|21|) and l|23|l we obtain the desired 
expression for the energy functional l|15|) in terms of the 
scattering data as 

= -h[r(k), {n n }} + -I 2 [r(k), {n n }} - I 3 [r(k)} 

7T 



K 



E(v„ \ M ! m ■ ■ h„ 

[jr ~ 1 ) Un + m- — logf 

n=l L 



A 



TO + K n 



(24) 



where 



oo 

tm-^ I '"^'r'^' U . (25) 
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We shall now extremize (|24() with respect to the scat- 
tering data, to obtain the self-consistent static fermion 
bags in the MGN model. Let us vary with respect to 
r{k) (with k > 0) first. As is evident from JTSJ, (t2*T)l and 
C3, 5£ re 9[a{x)]/5r{k) = F{k) r*(k)/(l - \r(k)\ 2 ), where 
F(k) is a calculable function, which does not vanish iden- 
tically. (For example, it can be shown that ImF(fc) = 
-(M/\ir)k/(m 2 + k 2 ).) Thus, r(k) = is the unique so- 
lution of the variational equation 5£ re9 [a(x)]/6r(k) = 0. 
Static extremal bags in the MGN model are reflection- 
less, as their counterparts in the GN model. 

Explicit formulas for reflectionless a(x) configurations 
with an arbitrary number K of p airs of bound states are 
displayed in Appendix B of 10]. In particular, the one 
which supports a single pair of bound states at energies 



±UJ b (K = ,/t 



the one originally discovered by 



DHN, is 



0~(x) = Tfl -(- K 



1 m + K 

tann kx log 

4 m — k 



tanh 



1 m + n 
- log 

4 m — k 



(26) 



We see that the formidable problem of finding the ex- 
tremal o~{x) configurations of the energy functional £ [a] 
(|15|l . is reduced to the simpler problem of extremizing 
an ordinary function £{uj u , c„) = £ [o~(x; co n , c n )\ with re- 
spect to the 2K parameters {c ni uj n } that determine the 
reflectionless background a(x). If we solve this ordinary 
extremum problem, we will be able to calculate the mass 
of the fermion bag. Let us write down this function ex- 
plicitly: 



£ re9 (Un) 



K 

E 

n=l 



2 2, /2Ae 

- + -log 

A 7r V m 



i | -T7 — 1 H — arctan — | ..c,, 

N 7T Kn 



M ( m — k,, 
— log — — 

A \m + K ri 



(27) 



The logarithmically divergent term in l|27|l should remind 
us that this equation is expressed in terms of the bare 
couplings A and M. As it turns out, the renormalization 
procedure of the effective potential for the constant con- 
densate, which we reviewed in the Introduction, applies 
also in the case of inhomogeneous background a(x): The 
bare coupling A in l|27(l is related to the renormalized one 
at the energy scale fi2 = m via an equation identical to 
©, namely, 



1 



1 



_ 1 2A 
A(A) A(to) 7r ^ m 



(28) 



(Due to the anisotropic cutoff implied in J3TJ, the cutoff 
scale in (|28f) is /«i = 2A rather than just A.) In addition, 
since the ratio ^4- is a renormalization group invariant 



all the way up to the cutoff scale, we can replace the 
coefficient of the last term in l|27|) by the common value 
of that invariant, namely, A ^ . With the help of these 
two relations, we obtain the renormalized form of (|27(l as 



K r 



A(m) 7T 



V n 2 LOft 

1 H — arctan — | uj n 

N 7T Kn 



m 



X(m] 



log 



m — k t 



m + K r 



(29) 



Next, note that the energy functional £ ren [a] evalu- 
ated at a reflectionless o~(x; uj n , c„), is independent of the 
normalization coefficients c„, that do affect the shape of 
a(x). The c„'s are thus flat directions of £ ren [o~] in the 
space of all reflectionless a(x) configurations. In fact, 
the c„'s (or more precisely, their logarithms) are collec- 
tive translational coordinates of the fermion bag o~(x) 
(see e.g., Appendix B in |l0()- One of these coordinates, 
corresponds, of course, to global translations of the bag 
as a whole. 

Finally, we are left with the task of extremizing l|29() 
with respect to the bound state energies u> n of the reflec- 
tionless background a(x). To this end, we follow DHN 
and parametrize these energies as u> n = mcos 9 n , with 
< 9 n < it/2. (Note that 8 n cannot attain the end-point 
values of its range: cu n = m plunges into the continuum, 
and u) n = is possible only if o~(x) is topologically non- 
trivial, which is not the case in the MGN model.) From 
t| 1 7JI we see that n n = m sin 9 n . In terms of the angular 
variables, we may write (|29|l as 



h«»)=»»e{(- 

71=1 ^ ^ 



A(m) 7T 



2" 

— I sm( 



V n 20 n 



COS 9 r , 



-7 log 



N 7T 

1 — sin 9 n 



1 + sin 9 n 



(30) 



where we have defined the renormalization group invari- 
ant 



f 
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M(/i) 



A(m) mA(/i) 
Finally, extremizing (|30|) with respect to 9 n yields 

ggren 



(31) 



d9 ri 



2m 



9 n v n 
~ ~ 2N 



7 tan ( 



sin 9 n = . 



(32) 



thus fixing 9 n as a function of the filling fraction 

< x n < 1 (33) 



N 



6 



according to 



h 7 tan n = — 

TT Z 



(34) 



The fact that the extremal value of 6 n is determined by 
the total number v n of particles and holes trapped in the 
bound states of the Dirac equation at ±u„, and not by 
the numbers of trapped particles and holes separately, is a 
manifestation of the underlying 0(2N) symmetry, which 
treats particles and holes symmetrically. As explained in 
(lif and in Appendix D of this fact indicates that 
this pair of bound states gives rise to an 0(2N) anti- 
symmetric tensor multiplet of rank v n of soliton states. 
The soliton as a whole is therefore the tensor product of 
all these antisymmetric tensor multiplets. The fermion 
number operator Nf is one of the generators of 0(2N). 
Its expectation value in the background of the extremal 
fermion bag in one of its 0(2N) states is (see Section 3 
of 0) 



K 



( N f) = Yl (p» - ■ 



(35) 



which is simply the sum over the individual valence 
fermion numbers 



N 



(«) 

f,val 



Pr, 



(36) 



of each of the states in the individual antisymmetric fac- 
tors. Evidently, for each of these antisymmetric repre- 
sentations 



< N {n) < v 



(37) 



in accordance with charge conjugation invariance. 

The left-hand side of (|34|l is a monotonically increas- 
ing function. Therefore, 134JI has a unique solution in 
the interval [0, tt/2\. This solution is evidently smaller 
than = the corresponding value of 8 n in the 
GN model for the same occupation number. Thus, the 
corresponding bound state energy u> n — mcos9 n in the 
MGN model is higher than its GN counterpart, and thus 
less bound. 

Substituting the extremal #'s from l|34l) in (|30|l we find 
that the mass Ai of our soliton (namely, N£ ren evaluated 
at the extremal point) is 



M({v„}) = Nm 



K 

E 

n=l 



2 . 
— sin 

TT 



9 n + 7 log 



1 + sin 6 r 
1 — sin r . 



(38) 

In the case of a single pair of bound states, K = 1, 13 2|) - 
(|34|l and l|38|) coincide with the corresponding results of 
variational calculations presented in |12l |13j , which were 
based on 1|26H as a trial configuration. In fact, it was 
realized in [13j that the trial configuration l|26|) is an exact 
solution of the extremum condition |f|jy = 0, provided 



(|34|l is used to fix k. This choice of trial configuration 
was very successful indeed! 

We should mention that renormalization of the energy 
functional l|15f) in the background of a generic reflection- 
less cr(x) and its extremization with respect to the n s 
can be carried out using an alternative method based on 
the diagonal resolvent of the Dirac operator ^j, which 
is basically a generalization of the calculations in 0, 0] 
to the case of an arbitrary number K of pairs of bound 
states. As it turns out, there are simple explicit formulas 
for the diagonal resolvents of the Dirac operator and of 
/i_ in reflectionless a(x) backgrounds, which make these 
computations possible, and lead to (|32|) - l|38|l fl9| . 



INVESTIGATING STABILITY OF EXTREMAL 
STATIC FERMION BAGS 

The extremal static soliton multiplets which we en- 
countered in the previous section, correspond, in the limit 
N — > oo, to exact eigenstates of the hamiltonian of the 
MGN model. However, at large but finite N, we expect 
some of these states to become unstable and thus to ac- 
quire small widths, similarly to the behavior of baryons in 
QCD with a large number of colors [2(J. The latter are 
also solitonic objects and are analogous to the "multi- 
quark" bound states of the MGN and GN models. In 
particular, Section 9 of [2(J offers a sketch of the l/N ex- 
pansion of two-dimensional QCD in the Coulomb gauge, 
both in the baryon and meson sectors, which is similar to 
the corresponding l/N expansion of the GN and MGN 
models, in the presence of fermion bags. (This expansion 
is based on the so-called bilocal condensate formalism, 
which was later developed in [tIEUE^.) 

Furthermore, we can imagine perturbing the MGN ac- 
tion by a small 0(2N) singlet perturbation (e.g., by 
adding to JIJ a term AS„ — e J d 2 xa n ), and ask which 
of the extremal fermion bags of the previous section are 
stable against such perturbations. (The perturbations 
Si and S2 correspond merely to a redefinition of the 
bare quantities M and <? 2 , and are thus not interest- 
ing. The higher perturbations, with n > 2, are non- 
renormalizablc. However, we could think of the resulting 
perturbed lagrangian as an effective one.) Under these 
circumstances, all possible decay channels of a given un- 
stable soliton multiplet must conserve energy, momentum 
and 0{2N) quantum numbers. 

It turns out that non-trivial results concerning stability 
may be established without getting into all the details of 
decomposing 0(2N) representations, by imposing a sim- 
ple necessary condition on the spectrum of the fermion 
number operator Nf in the multiplets involved in a given 
decay channel. This way of arguing (as described in de- 
tail below) has led, in the case of the GN model (see Sec- 
tion 4 of [lOj ) , to specification of all topologically-trivial 
stable fermion bags, consistent with the exact results of 
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[Tl| . Thus, it is reasonable to expect (and conjecture) 
that applying these stability considerations to the MGN 
model should lead to the correct list of fermion bags in 
this model which remain stable at finite N as well. Un- 
fortunately, exact results analogous to 01 an d valid for 
finite N, are not available for the MGN model, so this 
conjecture has to be verified by explicit calculation of 
1/7V corrections. We shall now make these stability con- 
siderations explicit. 

As we have learned so far, a given static soliton mul- 
tiplet in our model is a direct product of 0(2N) anti- 
symmetric tensors. The decay products of this soliton 
also correspond to a direct product of antisymmetric ten- 
sors. According to (|3*7|l . the spectrum of Nf in an an- 
tisymmetric tensor representation is symmetric, namely, 
_jymax < TV/ < iV™ a2: . When we compose two such rep- 
resentations Di,D 2 , the spectrum of Nf in the compos- 
ite representation D\ <g> D 2 , will obviously have the range 
-Nf ax {D 1 )-Nf ax {D 2 ) < N f (D 1 ®D 2 ) < N , f nax (D 1 ) + 
NJ lax (D 2 ). In particular, each of the possible eigenvalues 
in this range, will appear in at least one irreducible repre- 
sentation in the decomposition of Di <S> D 2 . More gener- 
ally, the spectrum of Nf{D\ ® D 2 ■ ■ ■ eg) D£) will have the 
range \Nf{D x ® ■■■^D L ) | < N^ ax (D 1 ) + - ■ ■ NJ iax (D L ). 

Consider now a decay process, in which a parent static 
soliton, which belongs to a (possibly reducible) represen- 
tation -Dparont, decays into a bunch of other solitons, such 
that the collection of all irreducible representations asso- 
ciated with the decay products is {D\, . . . , D^} (in which 
a given irreducible representation may occur more than 
once). By 0(2N) symmetry, the representation D pa rent 
must occur in the decomposition of D\ ® D 2 ■ ■ ■ (g> D^. 
Thus, according to the discussion in the previous para- 
graph, if this decay process is allowed, we must have 

Nf ax (L> parent ) < Nf ax {D x ) + ■■■ + N) nax (D L ) , (39) 

which is the necessary condition for 0(2N) symmetry 
we sought for. (Obviously, similar necessary conditions 
arise for the other N — 1 components of the highest weight 
vectors of the representations involved.) The decay pro- 
cess under consideration must respect energy conserva- 
tion, i.e., it must be exothermic. Thus, we supplement 
(|3"9"|) by the requirement 



M 



parent 



> 



products 



(40) 



on the masses Mi of the particles involved. 

For each of the static soliton multiplets discussed in the 
previous section, we will scan through all decay channels 
(into static solitons) and check which of these decay chan- 
nels are necessarily closed, simply by requiring that the 
two conditions (|39|) and H4()|l be mutually contradictory. 

In order to complete our argument, we must make the 
plausible assumption that for any time dependent stable 
soliton with given 0(2N) quantum numbers (should such 



a soliton exist), the static soliton with the same 0(2N) 
quantum numbers is lighter. In the GN this is in fact 
true for the known time dependent DHN breathers Q. 

For all the multiplets with a single pair of bound state, 
K = 1, we will find in this way that all decay channels 
are necessarily closed, thus establishing their stability. 
That these are stable multiplets is almost obvious to be- 
gin with - they are the lightest solitons, given their 0(2N) 
quantum numbers. We cannot establish in this way that 
all decay channels are necessarily closed for the higher 
solitons K > 1, and they are presumably unstable. 

Consider the function 

e{x) = -sm0(x) +7log , , < x < 1 , 

7r 1 — sm 0(x) 

(41) 

where 0{x) is a solution of pi^) . From and (|4*T)> we 

obtain that 



de(x) 
dx 



cos9(x) > 



d 2 e(x) it sin 9(x) 



dx 2 



2 1 + 7T7 sec 2 9(x] 



< 0. 



(42) 



Thus, e(x) is a monotonically increasing convex function 
in the range of interest, satisfying e(xi + x 2 ) < e(xi) + 
e(x 2 ). In terms of (|41|) . we may write the soliton mas 
simply as 



M(xx, - ■ ■ ,x n ) 



K 

Nm e(x r , 



(43) 



Now, we are ready to start the stability analysis. Con- 
sider the decaying parent soliton to be a static soli- 
ton with K pairs of bound states, corresponding to the 
direct product of K antisymmetric tensor representa- 
tions of ranks v\ , . . . , vk ■ The mass of this soliton is 
Ai(xi, . . . ,Xk) — ^ m ^2n=i e (^n), and according to 
(|37|l . the maximal fermion number eigenvalue occurring 
in this representation is NJ lax (D palcnt ) = J2n=i 

Following the strategy which we laid above, we shall 
now scan through all imaginable decay channels of this 
parent soliton (into final states of purely static solitons), 
and identify those channels which are necessarily closed. 
Thus, assume that the parent soliton under consideration 
decays into a configuration of lighter solitons, with quan- 
tum numbers of the direct product of L antisymmetric 
tensor representations v\ , . . . , vl- The way in which these 
L multiplets are arranged into extremal fermion bags is 
of no consequence to our discussion. Thus, we discuss all 
decay channels consistent with these quantum numbers 
in one sweep. 

The necessary conditions (|39|l and H4U|) for possible de- 
cay imply 



A" 



A' 



J2 in - ■ 



(44) 
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where all < x n , Xi < 1. The two pairs of boundary 
hypersurfaces in 144(1 are 

Si, Si : xi H h Xl = x\ H Y xr 

E 2 ,£ 2 : e(xi) H h e(x L ) = e(xi) H h e(5if ) (• 

where £i )2 are hypersurfaces in x-space and £1,2 are the 
corresponding hypersurfaces in i-space. 

Consider the behavior of (|43|) over the hyperplane 

S rct : xi H h Xk = T + a , (46) 

where < r < K is an integer and < a < 1. This hy- 
perplane corresponds, of course, to solitons with a fixed 
value of the maximal fermion number -/Vy" aa; (£'p aron t). 
We will now prove that M.(x\, ■ ■ ■ ,x n ) attains its abso- 
lute minimum on £ ra in the positive orthant, at the ver- 
tices of the intersection of E rct and the hypercube [0, 1]^, 
namely, the points 

X n ^ (f^nni ~1~ $nni "T" • ■ ■ ~\~ ^nnr)~^~^^nn r ^\ , — -!)•••; 

(47) 

with all possible choices of r + 1 coordinates ix, . . . , z r +i 
out of K . 

We prove this as follows (see Appendix F of [lOj): con- 
sider a sequence < x\ < ... < xk < 1, subjected 
to flin|l. Assume that for some i, < ij < Xi+i < 1. 
We will show that there exists another sequence of x's, 
with the same sum, but with a lower sum of the es. 
Thus, let 5 > be chosen such that Xi — 5 > and 
ajj+i + S < 1, i.e., < <5 < min{ii,l — ii+i}. Mod- 
ify the sequence under consideration by replacing Xi 
by Xi — 5 and Xi+\ by i^+i + 6, keeping the other 
K — 2 terms unaltered. The new sequence thus ob- 
tained has the same sum as the original sequence, and 
thus defines another point on E ra . We must show that 
D{8) — Aboriginal sequence) — A4(new sequence) > 0. 
Indeed, D(8) = e(xi) - e(xi - S) + e{x l+1 ) - e(x i+1 + S). 
Clearly, D(Q) — 0, and also D'(5) > in the relevant 
range of S. Thus, D(5) increases monotonically with 5, 
and reaches its maximum at 5 max = minjii, 1 — Xj+i}, 
where, depending on the initial condition at 5 = 0, cither 
Xi — <5 ma x = or Xi + i + <5 max = 1. Thus, the sequence 
of x's constrained to S rQ , which minimizes (|43(l . can- 
not have more than one element in the interior of [0, 1]. 
Thus, due to l|46|l . the absolute minimum is the sequence 
in which the r largest i's are 1, one x is a and the rest 
are zero, namely, the vertices (14711 . This is just the state- 
ment that the mass function ((43(1 , being the sum of the 
convex functions e(x), is convex inside the cube [0, 1]^. 

Therefore, a parent soliton corresponding to a point 
in the interior of the intersection of S rQ , and the hyper- 
cube [0, 1]^, can decay into a final state with quantum 
numbers corresponding to the points 1(47(1 . i.e., L = K 
and x n — Xn in ((44(1 . In fact, by continuity, such a par- 
ent soliton can decay also at least into the set of final 
states contained in small pockets above the vertices (14 711 , 



which correspond to L — K and x n — in' + e\ in 1(44(1 . 
where t\ is some small calculable quantity, or into final 
states corresponding to L > K in ((4411 . with Xi = x\ v ' 
for 1 < i < K and Xi = ti for K + 1 < i < L, and where 
£2 is another small calculable quantity. 

On the other hand, the parent soliton which corre- 
sponds to the vertices ((47(1 has no open channel to decay 
through. Hence it is potentially stable. Indeed, if it could 
decay through a channel corresponding to x\, . . . , Xl, 
then, from the requirement that these parameters sat- 
isfy ((44(1 . we would have 

L 

Xi > a + r 

i=l 

L 

< re(l)+e(a). (48) 

Define the hyperplane 

S ra : x% H h xl = a + r . (49) 

From the analysis above we know that the absolute min- 
imum of y]f—i € { x i) over the intersection of S ra and 
the hypercube [0, 1] L is re(l) + e(a). Thus, the points 
which satisfy the second inequality in ((48(1 are bounded 
by y'., Xi < a + r, in contradiction with the first in- 
equality in ((48(1 . This completes the proof that parent 
solitons which correspond to the vertices 1(47(1 are stable. 

Each such vertex represents a soliton which cannot de- 
cay through the channel under consideration, and is thus 
potentially stable. More precisely, all these vertices corre- 
spond to the same soliton, since the coordinates of these 
vertices are just permutations of each other, and thus 
all of them correspond to the same set of parameters, in 
which 

r of the x's are degenerate and equal to 1 

one of the x's is equal to a , and 

the remaining K — (r + 1) x's are null. (50) 

Does such a soliton exist? To answer this question 
let us recall a few basic facts: The parent soliton un- 
der discussion is topologically trivial. As such, it must 
bind fermions to be stabilized, and none of its bound 
state energies may vanish. Thus, all the ranks occur- 
ring in it must satisfy < v n < N. Finally, note, that 
due to the elementary fact, that the spectrum of the one 
dimensional Schrodinger operator /i_ cannot be degener- 
ate, all the w n 's must be different from each other, and 
so must be the v n s. Thus, the only physically realiz- 
able parent solitons, which are necessarily stable against 
the decay channel in question, correspond to r — and 
K = r + 1 = 1. These are, of course, the static solitons 
studied in [HEl 

A corollary of our analysis so far is that these K = 1 
solitons are stable against decaying into L f free fermions 
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or antifermions (i.e., Lf fundamental 0(2N) representa- 
tions) plus L — Lf solitons corresponding to higher an- 
tisymmetric tensor representations. In particular, it is 
stable against complete evaporation into free fermions. 
(Strictly speaking, this argument is valid only for values 
of Lf which are a finite fraction of TV, since our mass 
formula (|43|l is the leading order in the 1/N expansion, 
while removing a finite number of particles from the par- 
ent soliton is a perturbation of the order 1 /N relative to 
iffify .) Stability of the K = 1 solitons was studied in 
in detail. In particular, the regimes 7 << 1 and 7 >> 1 
were investigated. It was noted in |l2( that the mass of 
the most stable soliton, obtained as the filling fraction 
x — y 1, is nonanalytic around 7 = (with leading non- 
analytic behavior 7 log 7) . It was shown in ^3] that this 
nonanalyticity was related to the enhanced Z2 symme- 
try of the GN model at 7 = 0. The opposite regime, 
7 = M/Xm >> 1 may be attained by making the four- 
fermi interactions weak. In this regime the theory should 
be described in terms of quasi-free massive fermions of 
mass m. We thus expect that the binding energy of bags 
will tend to zero as 7 — * 00, which was indeed verified in 
d. 



APPENDIX: DERIVATION OF lt23l) 

Consider Schrodinger equation 

{-dl + V{x) - k 2 -m 2 ) V(x,k) = (51) 

with V(x) given by l|14|l. Due to the boundary conditions 
on <j{x) at spatial infinity, V(x) to 2 . Let <p(x, k) be 

its scattering solution 

t(k)4>(x,k) = t(k)e lkx +o(l) x ^+00 

t(k)(j)(x,k) = e lkx +r{k)e-' lkx +o(l) x -> -o<52) 

where t(k) and r(k) are, respectively, the transmission 
and reflection amplitudes. As a consequence of analytic- 
ity of these amplitudes, the transmission amplitude t(k) 
is completely determined on the real k axis by r{k) as 
141 



K 



t(k) = Vi-kWI 2 n 



IKl 



\l=l 



k — ini 



1 f log \l-\r(q)\ 2 } 

^ ' 2S PP - J q -k dq 



/ K 

n 

VJ=1 



IKl 



exp 



k — iki 



2ni 



log [1- |r(g)| s 
q — k — ie 



-dq 



(53) 



Evidently, the last expression can be extended to the 
upper-half complex fc-plane (as usual, e — y 0+), display- 
ing the bound states poles of t(k) in the upper half-plane 
explicitly. The function 



X 

^o(x) = c exp - J a{y)dy 



(54) 



solves (15 111 for uj = 0, i.e., for k = im. Since 
a(x) — to, it is not normalizable . Moreover, 
we can always choose the constant c such that 
y (x) — > e~ mx = e z(lm ^ x . Thus, with this choice of 

X—y-\-OC 

c, 



^o(x) = 4>{x 1 im) . 
Therefore, from Q52p. we obtain that 

(j)(L,im) 

t{im)e 



(55) 



2mL 



4>{-L,im) L^oo 

With the help of (|54|l we see that H56fl is equivalent to 



(56) 



{(j(x) — to) dx = — logi(im) . 



(57) 



Thus, finally, 



(<t(x) — to) dx = 
1 ? log [1- Kg)! 2 ] 



2ni 



K 



log 



im — q 



m — K r 



to + K Tl 



dq 



(58) 



from l|53|) . thus proving (|23|l . Due to the fact that on the 
real g-axis r(—q) — r*(q) and \t(q)\ 2 = 1 — \r(q)\ 2 < 1, 
we can rewrite (|58|l as 



(cr(x) — to) 



to f log |t(g)| 2 



A" 



7T / O 2 + TO 2 ' z — ' " V TO + Kr, 

" =1 



(59) 



JF thanks Michael Thies for useful correspondence. 
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